\(\int \frac {\sqrt {c+d x} \sqrt {e+f x} \sqrt {g+h x}}{a+b x} \, dx\) [43]

   Optimal result
   Rubi [A] (verified)
   Mathematica [C] (verified)
   Maple [A] (verified)
   Fricas [F(-1)]
   Sympy [F]
   Maxima [F]
   Giac [F]
   Mupad [F(-1)]

Optimal result

Integrand size = 35, antiderivative size = 570 \[ \int \frac {\sqrt {c+d x} \sqrt {e+f x} \sqrt {g+h x}}{a+b x} \, dx=\frac {2 \sqrt {c+d x} \sqrt {e+f x} \sqrt {g+h x}}{3 b}-\frac {2 \sqrt {-d e+c f} (3 a d f h-b (d f g+d e h+c f h)) \sqrt {\frac {d (e+f x)}{d e-c f}} \sqrt {g+h x} E\left (\arcsin \left (\frac {\sqrt {f} \sqrt {c+d x}}{\sqrt {-d e+c f}}\right )|\frac {(d e-c f) h}{f (d g-c h)}\right )}{3 b^2 d \sqrt {f} h \sqrt {e+f x} \sqrt {\frac {d (g+h x)}{d g-c h}}}+\frac {2 \sqrt {-d e+c f} \left (3 a^2 d f h^2-3 a b (d e+c f) h^2-b^2 (d g (f g-e h)-c h (f g+2 e h))\right ) \sqrt {\frac {d (e+f x)}{d e-c f}} \sqrt {\frac {d (g+h x)}{d g-c h}} \operatorname {EllipticF}\left (\arcsin \left (\frac {\sqrt {f} \sqrt {c+d x}}{\sqrt {-d e+c f}}\right ),\frac {(d e-c f) h}{f (d g-c h)}\right )}{3 b^3 d \sqrt {f} h \sqrt {e+f x} \sqrt {g+h x}}-\frac {2 (b e-a f) \sqrt {-d e+c f} (b g-a h) \sqrt {\frac {d (e+f x)}{d e-c f}} \sqrt {\frac {d (g+h x)}{d g-c h}} \operatorname {EllipticPi}\left (-\frac {b (d e-c f)}{(b c-a d) f},\arcsin \left (\frac {\sqrt {f} \sqrt {c+d x}}{\sqrt {-d e+c f}}\right ),\frac {(d e-c f) h}{f (d g-c h)}\right )}{b^3 \sqrt {f} \sqrt {e+f x} \sqrt {g+h x}} \]

[Out]

2/3*(d*x+c)^(1/2)*(f*x+e)^(1/2)*(h*x+g)^(1/2)/b-2/3*(3*a*d*f*h-b*(c*f*h+d*e*h+d*f*g))*EllipticE(f^(1/2)*(d*x+c
)^(1/2)/(c*f-d*e)^(1/2),((-c*f+d*e)*h/f/(-c*h+d*g))^(1/2))*(c*f-d*e)^(1/2)*(d*(f*x+e)/(-c*f+d*e))^(1/2)*(h*x+g
)^(1/2)/b^2/d/h/f^(1/2)/(f*x+e)^(1/2)/(d*(h*x+g)/(-c*h+d*g))^(1/2)+2/3*(3*a^2*d*f*h^2-3*a*b*(c*f+d*e)*h^2-b^2*
(d*g*(-e*h+f*g)-c*h*(2*e*h+f*g)))*EllipticF(f^(1/2)*(d*x+c)^(1/2)/(c*f-d*e)^(1/2),((-c*f+d*e)*h/f/(-c*h+d*g))^
(1/2))*(c*f-d*e)^(1/2)*(d*(f*x+e)/(-c*f+d*e))^(1/2)*(d*(h*x+g)/(-c*h+d*g))^(1/2)/b^3/d/h/f^(1/2)/(f*x+e)^(1/2)
/(h*x+g)^(1/2)-2*(-a*f+b*e)*(-a*h+b*g)*EllipticPi(f^(1/2)*(d*x+c)^(1/2)/(c*f-d*e)^(1/2),-b*(-c*f+d*e)/(-a*d+b*
c)/f,((-c*f+d*e)*h/f/(-c*h+d*g))^(1/2))*(c*f-d*e)^(1/2)*(d*(f*x+e)/(-c*f+d*e))^(1/2)*(d*(h*x+g)/(-c*h+d*g))^(1
/2)/b^3/f^(1/2)/(f*x+e)^(1/2)/(h*x+g)^(1/2)

Rubi [A] (verified)

Time = 0.87 (sec) , antiderivative size = 570, normalized size of antiderivative = 1.00, number of steps used = 12, number of rules used = 10, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.286, Rules used = {167, 1621, 175, 552, 551, 164, 115, 114, 122, 121} \[ \int \frac {\sqrt {c+d x} \sqrt {e+f x} \sqrt {g+h x}}{a+b x} \, dx=\frac {2 \sqrt {c f-d e} \sqrt {\frac {d (e+f x)}{d e-c f}} \sqrt {\frac {d (g+h x)}{d g-c h}} \left (3 a^2 d f h^2-3 a b h^2 (c f+d e)-\left (b^2 (d g (f g-e h)-c h (2 e h+f g))\right )\right ) \operatorname {EllipticF}\left (\arcsin \left (\frac {\sqrt {f} \sqrt {c+d x}}{\sqrt {c f-d e}}\right ),\frac {(d e-c f) h}{f (d g-c h)}\right )}{3 b^3 d \sqrt {f} h \sqrt {e+f x} \sqrt {g+h x}}-\frac {2 (b e-a f) (b g-a h) \sqrt {c f-d e} \sqrt {\frac {d (e+f x)}{d e-c f}} \sqrt {\frac {d (g+h x)}{d g-c h}} \operatorname {EllipticPi}\left (-\frac {b (d e-c f)}{(b c-a d) f},\arcsin \left (\frac {\sqrt {f} \sqrt {c+d x}}{\sqrt {c f-d e}}\right ),\frac {(d e-c f) h}{f (d g-c h)}\right )}{b^3 \sqrt {f} \sqrt {e+f x} \sqrt {g+h x}}-\frac {2 \sqrt {g+h x} \sqrt {c f-d e} \sqrt {\frac {d (e+f x)}{d e-c f}} (3 a d f h-b (c f h+d e h+d f g)) E\left (\arcsin \left (\frac {\sqrt {f} \sqrt {c+d x}}{\sqrt {c f-d e}}\right )|\frac {(d e-c f) h}{f (d g-c h)}\right )}{3 b^2 d \sqrt {f} h \sqrt {e+f x} \sqrt {\frac {d (g+h x)}{d g-c h}}}+\frac {2 \sqrt {c+d x} \sqrt {e+f x} \sqrt {g+h x}}{3 b} \]

[In]

Int[(Sqrt[c + d*x]*Sqrt[e + f*x]*Sqrt[g + h*x])/(a + b*x),x]

[Out]

(2*Sqrt[c + d*x]*Sqrt[e + f*x]*Sqrt[g + h*x])/(3*b) - (2*Sqrt[-(d*e) + c*f]*(3*a*d*f*h - b*(d*f*g + d*e*h + c*
f*h))*Sqrt[(d*(e + f*x))/(d*e - c*f)]*Sqrt[g + h*x]*EllipticE[ArcSin[(Sqrt[f]*Sqrt[c + d*x])/Sqrt[-(d*e) + c*f
]], ((d*e - c*f)*h)/(f*(d*g - c*h))])/(3*b^2*d*Sqrt[f]*h*Sqrt[e + f*x]*Sqrt[(d*(g + h*x))/(d*g - c*h)]) + (2*S
qrt[-(d*e) + c*f]*(3*a^2*d*f*h^2 - 3*a*b*(d*e + c*f)*h^2 - b^2*(d*g*(f*g - e*h) - c*h*(f*g + 2*e*h)))*Sqrt[(d*
(e + f*x))/(d*e - c*f)]*Sqrt[(d*(g + h*x))/(d*g - c*h)]*EllipticF[ArcSin[(Sqrt[f]*Sqrt[c + d*x])/Sqrt[-(d*e) +
 c*f]], ((d*e - c*f)*h)/(f*(d*g - c*h))])/(3*b^3*d*Sqrt[f]*h*Sqrt[e + f*x]*Sqrt[g + h*x]) - (2*(b*e - a*f)*Sqr
t[-(d*e) + c*f]*(b*g - a*h)*Sqrt[(d*(e + f*x))/(d*e - c*f)]*Sqrt[(d*(g + h*x))/(d*g - c*h)]*EllipticPi[-((b*(d
*e - c*f))/((b*c - a*d)*f)), ArcSin[(Sqrt[f]*Sqrt[c + d*x])/Sqrt[-(d*e) + c*f]], ((d*e - c*f)*h)/(f*(d*g - c*h
))])/(b^3*Sqrt[f]*Sqrt[e + f*x]*Sqrt[g + h*x])

Rule 114

Int[Sqrt[(e_.) + (f_.)*(x_)]/(Sqrt[(a_) + (b_.)*(x_)]*Sqrt[(c_) + (d_.)*(x_)]), x_Symbol] :> Simp[(2/b)*Rt[-(b
*e - a*f)/d, 2]*EllipticE[ArcSin[Sqrt[a + b*x]/Rt[-(b*c - a*d)/d, 2]], f*((b*c - a*d)/(d*(b*e - a*f)))], x] /;
 FreeQ[{a, b, c, d, e, f}, x] && GtQ[b/(b*c - a*d), 0] && GtQ[b/(b*e - a*f), 0] &&  !LtQ[-(b*c - a*d)/d, 0] &&
  !(SimplerQ[c + d*x, a + b*x] && GtQ[-d/(b*c - a*d), 0] && GtQ[d/(d*e - c*f), 0] &&  !LtQ[(b*c - a*d)/b, 0])

Rule 115

Int[Sqrt[(e_.) + (f_.)*(x_)]/(Sqrt[(a_) + (b_.)*(x_)]*Sqrt[(c_) + (d_.)*(x_)]), x_Symbol] :> Dist[Sqrt[e + f*x
]*(Sqrt[b*((c + d*x)/(b*c - a*d))]/(Sqrt[c + d*x]*Sqrt[b*((e + f*x)/(b*e - a*f))])), Int[Sqrt[b*(e/(b*e - a*f)
) + b*f*(x/(b*e - a*f))]/(Sqrt[a + b*x]*Sqrt[b*(c/(b*c - a*d)) + b*d*(x/(b*c - a*d))]), x], x] /; FreeQ[{a, b,
 c, d, e, f}, x] &&  !(GtQ[b/(b*c - a*d), 0] && GtQ[b/(b*e - a*f), 0]) &&  !LtQ[-(b*c - a*d)/d, 0]

Rule 121

Int[1/(Sqrt[(a_) + (b_.)*(x_)]*Sqrt[(c_) + (d_.)*(x_)]*Sqrt[(e_) + (f_.)*(x_)]), x_Symbol] :> Simp[2*(Rt[-b/d,
 2]/(b*Sqrt[(b*e - a*f)/b]))*EllipticF[ArcSin[Sqrt[a + b*x]/(Rt[-b/d, 2]*Sqrt[(b*c - a*d)/b])], f*((b*c - a*d)
/(d*(b*e - a*f)))], x] /; FreeQ[{a, b, c, d, e, f}, x] && GtQ[b/(b*c - a*d), 0] && GtQ[b/(b*e - a*f), 0] && Si
mplerQ[a + b*x, c + d*x] && SimplerQ[a + b*x, e + f*x] && (PosQ[-(b*c - a*d)/d] || NegQ[-(b*e - a*f)/f])

Rule 122

Int[1/(Sqrt[(a_) + (b_.)*(x_)]*Sqrt[(c_) + (d_.)*(x_)]*Sqrt[(e_) + (f_.)*(x_)]), x_Symbol] :> Dist[Sqrt[b*((c
+ d*x)/(b*c - a*d))]/Sqrt[c + d*x], Int[1/(Sqrt[a + b*x]*Sqrt[b*(c/(b*c - a*d)) + b*d*(x/(b*c - a*d))]*Sqrt[e
+ f*x]), x], x] /; FreeQ[{a, b, c, d, e, f}, x] &&  !GtQ[(b*c - a*d)/b, 0] && SimplerQ[a + b*x, c + d*x] && Si
mplerQ[a + b*x, e + f*x]

Rule 164

Int[((g_.) + (h_.)*(x_))/(Sqrt[(a_.) + (b_.)*(x_)]*Sqrt[(c_) + (d_.)*(x_)]*Sqrt[(e_) + (f_.)*(x_)]), x_Symbol]
 :> Dist[h/f, Int[Sqrt[e + f*x]/(Sqrt[a + b*x]*Sqrt[c + d*x]), x], x] + Dist[(f*g - e*h)/f, Int[1/(Sqrt[a + b*
x]*Sqrt[c + d*x]*Sqrt[e + f*x]), x], x] /; FreeQ[{a, b, c, d, e, f, g, h}, x] && SimplerQ[a + b*x, e + f*x] &&
 SimplerQ[c + d*x, e + f*x]

Rule 167

Int[((a_.) + (b_.)*(x_))^(m_)*Sqrt[(c_.) + (d_.)*(x_)]*Sqrt[(e_.) + (f_.)*(x_)]*Sqrt[(g_.) + (h_.)*(x_)], x_Sy
mbol] :> Simp[2*(a + b*x)^(m + 1)*Sqrt[c + d*x]*Sqrt[e + f*x]*(Sqrt[g + h*x]/(b*(2*m + 5))), x] + Dist[1/(b*(2
*m + 5)), Int[((a + b*x)^m/(Sqrt[c + d*x]*Sqrt[e + f*x]*Sqrt[g + h*x]))*Simp[3*b*c*e*g - a*(d*e*g + c*f*g + c*
e*h) + 2*(b*(d*e*g + c*f*g + c*e*h) - a*(d*f*g + d*e*h + c*f*h))*x - (3*a*d*f*h - b*(d*f*g + d*e*h + c*f*h))*x
^2, x], x], x] /; FreeQ[{a, b, c, d, e, f, g, h, m}, x] && IntegerQ[2*m] &&  !LtQ[m, -1]

Rule 175

Int[1/(((a_.) + (b_.)*(x_))*Sqrt[(c_.) + (d_.)*(x_)]*Sqrt[(e_.) + (f_.)*(x_)]*Sqrt[(g_.) + (h_.)*(x_)]), x_Sym
bol] :> Dist[-2, Subst[Int[1/(Simp[b*c - a*d - b*x^2, x]*Sqrt[Simp[(d*e - c*f)/d + f*(x^2/d), x]]*Sqrt[Simp[(d
*g - c*h)/d + h*(x^2/d), x]]), x], x, Sqrt[c + d*x]], x] /; FreeQ[{a, b, c, d, e, f, g, h}, x] &&  !SimplerQ[e
 + f*x, c + d*x] &&  !SimplerQ[g + h*x, c + d*x]

Rule 551

Int[1/(((a_) + (b_.)*(x_)^2)*Sqrt[(c_) + (d_.)*(x_)^2]*Sqrt[(e_) + (f_.)*(x_)^2]), x_Symbol] :> Simp[(1/(a*Sqr
t[c]*Sqrt[e]*Rt[-d/c, 2]))*EllipticPi[b*(c/(a*d)), ArcSin[Rt[-d/c, 2]*x], c*(f/(d*e))], x] /; FreeQ[{a, b, c,
d, e, f}, x] &&  !GtQ[d/c, 0] && GtQ[c, 0] && GtQ[e, 0] &&  !( !GtQ[f/e, 0] && SimplerSqrtQ[-f/e, -d/c])

Rule 552

Int[1/(((a_) + (b_.)*(x_)^2)*Sqrt[(c_) + (d_.)*(x_)^2]*Sqrt[(e_) + (f_.)*(x_)^2]), x_Symbol] :> Dist[Sqrt[1 +
(d/c)*x^2]/Sqrt[c + d*x^2], Int[1/((a + b*x^2)*Sqrt[1 + (d/c)*x^2]*Sqrt[e + f*x^2]), x], x] /; FreeQ[{a, b, c,
 d, e, f}, x] &&  !GtQ[c, 0]

Rule 1621

Int[(Px_)*((a_.) + (b_.)*(x_))^(m_.)*((c_.) + (d_.)*(x_))^(n_.)*((e_.) + (f_.)*(x_))^(p_.)*((g_.) + (h_.)*(x_)
)^(q_.), x_Symbol] :> Dist[PolynomialRemainder[Px, a + b*x, x], Int[(a + b*x)^m*(c + d*x)^n*(e + f*x)^p*(g + h
*x)^q, x], x] + Int[PolynomialQuotient[Px, a + b*x, x]*(a + b*x)^(m + 1)*(c + d*x)^n*(e + f*x)^p*(g + h*x)^q,
x] /; FreeQ[{a, b, c, d, e, f, g, h, m, n, p, q}, x] && PolyQ[Px, x] && EqQ[m, -1]

Rubi steps \begin{align*} \text {integral}& = \frac {2 \sqrt {c+d x} \sqrt {e+f x} \sqrt {g+h x}}{3 b}+\frac {\int \frac {3 b c e g-a (d e g+c f g+c e h)+2 (b (d e g+c f g+c e h)-a (d f g+d e h+c f h)) x-(3 a d f h-b (d f g+d e h+c f h)) x^2}{(a+b x) \sqrt {c+d x} \sqrt {e+f x} \sqrt {g+h x}} \, dx}{3 b} \\ & = \frac {2 \sqrt {c+d x} \sqrt {e+f x} \sqrt {g+h x}}{3 b}+\frac {\int \frac {2 d e g+2 c f g-\frac {3 a d f g}{b}+2 c e h-\frac {3 a d e h}{b}-\frac {3 a c f h}{b}+\frac {3 a^2 d f h}{b^2}+\left (d f g+d e h+c f h-\frac {3 a d f h}{b}\right ) x}{\sqrt {c+d x} \sqrt {e+f x} \sqrt {g+h x}} \, dx}{3 b}+\frac {((b c-a d) (b e-a f) (b g-a h)) \int \frac {1}{(a+b x) \sqrt {c+d x} \sqrt {e+f x} \sqrt {g+h x}} \, dx}{b^3} \\ & = \frac {2 \sqrt {c+d x} \sqrt {e+f x} \sqrt {g+h x}}{3 b}-\frac {(2 (b c-a d) (b e-a f) (b g-a h)) \text {Subst}\left (\int \frac {1}{\left (b c-a d-b x^2\right ) \sqrt {e-\frac {c f}{d}+\frac {f x^2}{d}} \sqrt {g-\frac {c h}{d}+\frac {h x^2}{d}}} \, dx,x,\sqrt {c+d x}\right )}{b^3}+\frac {\left (d f g+d e h+c f h-\frac {3 a d f h}{b}\right ) \int \frac {\sqrt {g+h x}}{\sqrt {c+d x} \sqrt {e+f x}} \, dx}{3 b h}+\frac {\left (h \left (2 d e g+2 c f g-\frac {3 a d f g}{b}+2 c e h-\frac {3 a d e h}{b}-\frac {3 a c f h}{b}+\frac {3 a^2 d f h}{b^2}\right )-g \left (d f g+d e h+c f h-\frac {3 a d f h}{b}\right )\right ) \int \frac {1}{\sqrt {c+d x} \sqrt {e+f x} \sqrt {g+h x}} \, dx}{3 b h} \\ & = \frac {2 \sqrt {c+d x} \sqrt {e+f x} \sqrt {g+h x}}{3 b}-\frac {\left (2 (b c-a d) (b e-a f) (b g-a h) \sqrt {\frac {d (e+f x)}{d e-c f}}\right ) \text {Subst}\left (\int \frac {1}{\left (b c-a d-b x^2\right ) \sqrt {1+\frac {f x^2}{d \left (e-\frac {c f}{d}\right )}} \sqrt {g-\frac {c h}{d}+\frac {h x^2}{d}}} \, dx,x,\sqrt {c+d x}\right )}{b^3 \sqrt {e+f x}}+\frac {\left (\left (h \left (2 d e g+2 c f g-\frac {3 a d f g}{b}+2 c e h-\frac {3 a d e h}{b}-\frac {3 a c f h}{b}+\frac {3 a^2 d f h}{b^2}\right )-g \left (d f g+d e h+c f h-\frac {3 a d f h}{b}\right )\right ) \sqrt {\frac {d (e+f x)}{d e-c f}}\right ) \int \frac {1}{\sqrt {c+d x} \sqrt {\frac {d e}{d e-c f}+\frac {d f x}{d e-c f}} \sqrt {g+h x}} \, dx}{3 b h \sqrt {e+f x}}+\frac {\left (\left (d f g+d e h+c f h-\frac {3 a d f h}{b}\right ) \sqrt {\frac {d (e+f x)}{d e-c f}} \sqrt {g+h x}\right ) \int \frac {\sqrt {\frac {d g}{d g-c h}+\frac {d h x}{d g-c h}}}{\sqrt {c+d x} \sqrt {\frac {d e}{d e-c f}+\frac {d f x}{d e-c f}}} \, dx}{3 b h \sqrt {e+f x} \sqrt {\frac {d (g+h x)}{d g-c h}}} \\ & = \frac {2 \sqrt {c+d x} \sqrt {e+f x} \sqrt {g+h x}}{3 b}-\frac {2 \sqrt {-d e+c f} (3 a d f h-b (d f g+d e h+c f h)) \sqrt {\frac {d (e+f x)}{d e-c f}} \sqrt {g+h x} E\left (\sin ^{-1}\left (\frac {\sqrt {f} \sqrt {c+d x}}{\sqrt {-d e+c f}}\right )|\frac {(d e-c f) h}{f (d g-c h)}\right )}{3 b^2 d \sqrt {f} h \sqrt {e+f x} \sqrt {\frac {d (g+h x)}{d g-c h}}}-\frac {\left (2 (b c-a d) (b e-a f) (b g-a h) \sqrt {\frac {d (e+f x)}{d e-c f}} \sqrt {\frac {d (g+h x)}{d g-c h}}\right ) \text {Subst}\left (\int \frac {1}{\left (b c-a d-b x^2\right ) \sqrt {1+\frac {f x^2}{d \left (e-\frac {c f}{d}\right )}} \sqrt {1+\frac {h x^2}{d \left (g-\frac {c h}{d}\right )}}} \, dx,x,\sqrt {c+d x}\right )}{b^3 \sqrt {e+f x} \sqrt {g+h x}}+\frac {\left (\left (h \left (2 d e g+2 c f g-\frac {3 a d f g}{b}+2 c e h-\frac {3 a d e h}{b}-\frac {3 a c f h}{b}+\frac {3 a^2 d f h}{b^2}\right )-g \left (d f g+d e h+c f h-\frac {3 a d f h}{b}\right )\right ) \sqrt {\frac {d (e+f x)}{d e-c f}} \sqrt {\frac {d (g+h x)}{d g-c h}}\right ) \int \frac {1}{\sqrt {c+d x} \sqrt {\frac {d e}{d e-c f}+\frac {d f x}{d e-c f}} \sqrt {\frac {d g}{d g-c h}+\frac {d h x}{d g-c h}}} \, dx}{3 b h \sqrt {e+f x} \sqrt {g+h x}} \\ & = \frac {2 \sqrt {c+d x} \sqrt {e+f x} \sqrt {g+h x}}{3 b}-\frac {2 \sqrt {-d e+c f} (3 a d f h-b (d f g+d e h+c f h)) \sqrt {\frac {d (e+f x)}{d e-c f}} \sqrt {g+h x} E\left (\sin ^{-1}\left (\frac {\sqrt {f} \sqrt {c+d x}}{\sqrt {-d e+c f}}\right )|\frac {(d e-c f) h}{f (d g-c h)}\right )}{3 b^2 d \sqrt {f} h \sqrt {e+f x} \sqrt {\frac {d (g+h x)}{d g-c h}}}+\frac {2 \sqrt {-d e+c f} \left (3 a^2 d f h^2-3 a b (d e+c f) h^2-b^2 (d g (f g-e h)-c h (f g+2 e h))\right ) \sqrt {\frac {d (e+f x)}{d e-c f}} \sqrt {\frac {d (g+h x)}{d g-c h}} F\left (\sin ^{-1}\left (\frac {\sqrt {f} \sqrt {c+d x}}{\sqrt {-d e+c f}}\right )|\frac {(d e-c f) h}{f (d g-c h)}\right )}{3 b^3 d \sqrt {f} h \sqrt {e+f x} \sqrt {g+h x}}-\frac {2 (b e-a f) \sqrt {-d e+c f} (b g-a h) \sqrt {\frac {d (e+f x)}{d e-c f}} \sqrt {\frac {d (g+h x)}{d g-c h}} \Pi \left (-\frac {b (d e-c f)}{(b c-a d) f};\sin ^{-1}\left (\frac {\sqrt {f} \sqrt {c+d x}}{\sqrt {-d e+c f}}\right )|\frac {(d e-c f) h}{f (d g-c h)}\right )}{b^3 \sqrt {f} \sqrt {e+f x} \sqrt {g+h x}} \\ \end{align*}

Mathematica [C] (verified)

Result contains complex when optimal does not.

Time = 30.02 (sec) , antiderivative size = 1254, normalized size of antiderivative = 2.20 \[ \int \frac {\sqrt {c+d x} \sqrt {e+f x} \sqrt {g+h x}}{a+b x} \, dx=\frac {2 \sqrt {c+d x} \left (3 b^2 e g-3 a b f g+\frac {b^2 f g^2}{h}-3 a b e h+\frac {b^2 e^2 h}{f}-\frac {b^2 c^2 f h}{d^2}+\frac {3 a b c f h}{d}+2 b^2 f g x+2 b^2 e h x-3 a b f h x+\frac {b^2 c f h x}{d}+b^2 f h x^2-\frac {b^2 c e g}{c+d x}-\frac {3 a b d e g}{c+d x}+\frac {3 a b c f g}{c+d x}+\frac {3 a b c e h}{c+d x}+\frac {b^2 c^3 f h}{d^2 (c+d x)}-\frac {3 a b c^2 f h}{d (c+d x)}+\frac {b^2 d e^2 g}{c f+d f x}-\frac {b^2 c e^2 h}{c f+d f x}+\frac {b^2 d e g^2}{c h+d h x}-\frac {b^2 c f g^2}{c h+d h x}-\frac {i b \sqrt {-c+\frac {d e}{f}} (3 a d f h-b (d f g+d e h+c f h)) \sqrt {c+d x} \sqrt {\frac {d (e+f x)}{f (c+d x)}} \sqrt {\frac {d (g+h x)}{h (c+d x)}} E\left (i \text {arcsinh}\left (\frac {\sqrt {-c+\frac {d e}{f}}}{\sqrt {c+d x}}\right )|\frac {d f g-c f h}{d e h-c f h}\right )}{d^2}+\frac {i b \sqrt {-c+\frac {d e}{f}} (-2 b f g-b e h+3 a f h) \sqrt {c+d x} \sqrt {\frac {d (e+f x)}{f (c+d x)}} \sqrt {\frac {d (g+h x)}{h (c+d x)}} \operatorname {EllipticF}\left (i \text {arcsinh}\left (\frac {\sqrt {-c+\frac {d e}{f}}}{\sqrt {c+d x}}\right ),\frac {d f g-c f h}{d e h-c f h}\right )}{d}+\frac {3 i b^2 e \sqrt {-c+\frac {d e}{f}} f g \sqrt {c+d x} \sqrt {\frac {d (e+f x)}{f (c+d x)}} \sqrt {\frac {d (g+h x)}{h (c+d x)}} \operatorname {EllipticPi}\left (-\frac {b c f-a d f}{b d e-b c f},i \text {arcsinh}\left (\frac {\sqrt {-c+\frac {d e}{f}}}{\sqrt {c+d x}}\right ),\frac {d f g-c f h}{d e h-c f h}\right )}{d e-c f}+\frac {3 i a b \sqrt {-c+\frac {d e}{f}} f^2 g \sqrt {c+d x} \sqrt {\frac {d (e+f x)}{f (c+d x)}} \sqrt {\frac {d (g+h x)}{h (c+d x)}} \operatorname {EllipticPi}\left (-\frac {b c f-a d f}{b d e-b c f},i \text {arcsinh}\left (\frac {\sqrt {-c+\frac {d e}{f}}}{\sqrt {c+d x}}\right ),\frac {d f g-c f h}{d e h-c f h}\right )}{-d e+c f}+\frac {3 i a^2 \sqrt {-c+\frac {d e}{f}} f^2 h \sqrt {c+d x} \sqrt {\frac {d (e+f x)}{f (c+d x)}} \sqrt {\frac {d (g+h x)}{h (c+d x)}} \operatorname {EllipticPi}\left (-\frac {b c f-a d f}{b d e-b c f},i \text {arcsinh}\left (\frac {\sqrt {-c+\frac {d e}{f}}}{\sqrt {c+d x}}\right ),\frac {d f g-c f h}{d e h-c f h}\right )}{d e-c f}+\frac {3 i a b e \sqrt {-c+\frac {d e}{f}} f h \sqrt {c+d x} \sqrt {\frac {d (e+f x)}{f (c+d x)}} \sqrt {\frac {d (g+h x)}{h (c+d x)}} \operatorname {EllipticPi}\left (-\frac {b c f-a d f}{b d e-b c f},i \text {arcsinh}\left (\frac {\sqrt {-c+\frac {d e}{f}}}{\sqrt {c+d x}}\right ),\frac {d f g-c f h}{d e h-c f h}\right )}{-d e+c f}\right )}{3 b^3 \sqrt {e+f x} \sqrt {g+h x}} \]

[In]

Integrate[(Sqrt[c + d*x]*Sqrt[e + f*x]*Sqrt[g + h*x])/(a + b*x),x]

[Out]

(2*Sqrt[c + d*x]*(3*b^2*e*g - 3*a*b*f*g + (b^2*f*g^2)/h - 3*a*b*e*h + (b^2*e^2*h)/f - (b^2*c^2*f*h)/d^2 + (3*a
*b*c*f*h)/d + 2*b^2*f*g*x + 2*b^2*e*h*x - 3*a*b*f*h*x + (b^2*c*f*h*x)/d + b^2*f*h*x^2 - (b^2*c*e*g)/(c + d*x)
- (3*a*b*d*e*g)/(c + d*x) + (3*a*b*c*f*g)/(c + d*x) + (3*a*b*c*e*h)/(c + d*x) + (b^2*c^3*f*h)/(d^2*(c + d*x))
- (3*a*b*c^2*f*h)/(d*(c + d*x)) + (b^2*d*e^2*g)/(c*f + d*f*x) - (b^2*c*e^2*h)/(c*f + d*f*x) + (b^2*d*e*g^2)/(c
*h + d*h*x) - (b^2*c*f*g^2)/(c*h + d*h*x) - (I*b*Sqrt[-c + (d*e)/f]*(3*a*d*f*h - b*(d*f*g + d*e*h + c*f*h))*Sq
rt[c + d*x]*Sqrt[(d*(e + f*x))/(f*(c + d*x))]*Sqrt[(d*(g + h*x))/(h*(c + d*x))]*EllipticE[I*ArcSinh[Sqrt[-c +
(d*e)/f]/Sqrt[c + d*x]], (d*f*g - c*f*h)/(d*e*h - c*f*h)])/d^2 + (I*b*Sqrt[-c + (d*e)/f]*(-2*b*f*g - b*e*h + 3
*a*f*h)*Sqrt[c + d*x]*Sqrt[(d*(e + f*x))/(f*(c + d*x))]*Sqrt[(d*(g + h*x))/(h*(c + d*x))]*EllipticF[I*ArcSinh[
Sqrt[-c + (d*e)/f]/Sqrt[c + d*x]], (d*f*g - c*f*h)/(d*e*h - c*f*h)])/d + ((3*I)*b^2*e*Sqrt[-c + (d*e)/f]*f*g*S
qrt[c + d*x]*Sqrt[(d*(e + f*x))/(f*(c + d*x))]*Sqrt[(d*(g + h*x))/(h*(c + d*x))]*EllipticPi[-((b*c*f - a*d*f)/
(b*d*e - b*c*f)), I*ArcSinh[Sqrt[-c + (d*e)/f]/Sqrt[c + d*x]], (d*f*g - c*f*h)/(d*e*h - c*f*h)])/(d*e - c*f) +
 ((3*I)*a*b*Sqrt[-c + (d*e)/f]*f^2*g*Sqrt[c + d*x]*Sqrt[(d*(e + f*x))/(f*(c + d*x))]*Sqrt[(d*(g + h*x))/(h*(c
+ d*x))]*EllipticPi[-((b*c*f - a*d*f)/(b*d*e - b*c*f)), I*ArcSinh[Sqrt[-c + (d*e)/f]/Sqrt[c + d*x]], (d*f*g -
c*f*h)/(d*e*h - c*f*h)])/(-(d*e) + c*f) + ((3*I)*a^2*Sqrt[-c + (d*e)/f]*f^2*h*Sqrt[c + d*x]*Sqrt[(d*(e + f*x))
/(f*(c + d*x))]*Sqrt[(d*(g + h*x))/(h*(c + d*x))]*EllipticPi[-((b*c*f - a*d*f)/(b*d*e - b*c*f)), I*ArcSinh[Sqr
t[-c + (d*e)/f]/Sqrt[c + d*x]], (d*f*g - c*f*h)/(d*e*h - c*f*h)])/(d*e - c*f) + ((3*I)*a*b*e*Sqrt[-c + (d*e)/f
]*f*h*Sqrt[c + d*x]*Sqrt[(d*(e + f*x))/(f*(c + d*x))]*Sqrt[(d*(g + h*x))/(h*(c + d*x))]*EllipticPi[-((b*c*f -
a*d*f)/(b*d*e - b*c*f)), I*ArcSinh[Sqrt[-c + (d*e)/f]/Sqrt[c + d*x]], (d*f*g - c*f*h)/(d*e*h - c*f*h)])/(-(d*e
) + c*f)))/(3*b^3*Sqrt[e + f*x]*Sqrt[g + h*x])

Maple [A] (verified)

Time = 1.64 (sec) , antiderivative size = 976, normalized size of antiderivative = 1.71

method result size
elliptic \(\frac {\sqrt {\left (d x +c \right ) \left (f x +e \right ) \left (h x +g \right )}\, \left (\frac {2 \sqrt {d f h \,x^{3}+c f h \,x^{2}+d e h \,x^{2}+d f g \,x^{2}+c e h x +c f g x +d e g x +c e g}}{3 b}+\frac {2 \left (\frac {a^{2} d f h -a b c f h -a b d e h -a b d f g +b^{2} c e h +b^{2} c f g +b^{2} d e g}{b^{3}}-\frac {2 \left (\frac {1}{2} c e h +\frac {1}{2} c f g +\frac {1}{2} d e g \right )}{3 b}\right ) \left (\frac {g}{h}-\frac {e}{f}\right ) \sqrt {\frac {x +\frac {g}{h}}{\frac {g}{h}-\frac {e}{f}}}\, \sqrt {\frac {x +\frac {c}{d}}{-\frac {g}{h}+\frac {c}{d}}}\, \sqrt {\frac {x +\frac {e}{f}}{-\frac {g}{h}+\frac {e}{f}}}\, F\left (\sqrt {\frac {x +\frac {g}{h}}{\frac {g}{h}-\frac {e}{f}}}, \sqrt {\frac {-\frac {g}{h}+\frac {e}{f}}{-\frac {g}{h}+\frac {c}{d}}}\right )}{\sqrt {d f h \,x^{3}+c f h \,x^{2}+d e h \,x^{2}+d f g \,x^{2}+c e h x +c f g x +d e g x +c e g}}+\frac {2 \left (-\frac {a d f h -b c f h -b d e h -b d f g}{b^{2}}-\frac {2 \left (c f h +d e h +d f g \right )}{3 b}\right ) \left (\frac {g}{h}-\frac {e}{f}\right ) \sqrt {\frac {x +\frac {g}{h}}{\frac {g}{h}-\frac {e}{f}}}\, \sqrt {\frac {x +\frac {c}{d}}{-\frac {g}{h}+\frac {c}{d}}}\, \sqrt {\frac {x +\frac {e}{f}}{-\frac {g}{h}+\frac {e}{f}}}\, \left (\left (-\frac {g}{h}+\frac {c}{d}\right ) E\left (\sqrt {\frac {x +\frac {g}{h}}{\frac {g}{h}-\frac {e}{f}}}, \sqrt {\frac {-\frac {g}{h}+\frac {e}{f}}{-\frac {g}{h}+\frac {c}{d}}}\right )-\frac {c F\left (\sqrt {\frac {x +\frac {g}{h}}{\frac {g}{h}-\frac {e}{f}}}, \sqrt {\frac {-\frac {g}{h}+\frac {e}{f}}{-\frac {g}{h}+\frac {c}{d}}}\right )}{d}\right )}{\sqrt {d f h \,x^{3}+c f h \,x^{2}+d e h \,x^{2}+d f g \,x^{2}+c e h x +c f g x +d e g x +c e g}}-\frac {2 \left (a^{3} d f h -a^{2} b c f h -a^{2} b d e h -a^{2} b d f g +a \,b^{2} c e h +a \,b^{2} c f g +a \,b^{2} d e g -b^{3} c e g \right ) \left (\frac {g}{h}-\frac {e}{f}\right ) \sqrt {\frac {x +\frac {g}{h}}{\frac {g}{h}-\frac {e}{f}}}\, \sqrt {\frac {x +\frac {c}{d}}{-\frac {g}{h}+\frac {c}{d}}}\, \sqrt {\frac {x +\frac {e}{f}}{-\frac {g}{h}+\frac {e}{f}}}\, \Pi \left (\sqrt {\frac {x +\frac {g}{h}}{\frac {g}{h}-\frac {e}{f}}}, \frac {-\frac {g}{h}+\frac {e}{f}}{-\frac {g}{h}+\frac {a}{b}}, \sqrt {\frac {-\frac {g}{h}+\frac {e}{f}}{-\frac {g}{h}+\frac {c}{d}}}\right )}{b^{4} \sqrt {d f h \,x^{3}+c f h \,x^{2}+d e h \,x^{2}+d f g \,x^{2}+c e h x +c f g x +d e g x +c e g}\, \left (-\frac {g}{h}+\frac {a}{b}\right )}\right )}{\sqrt {d x +c}\, \sqrt {f x +e}\, \sqrt {h x +g}}\) \(976\)
default \(\text {Expression too large to display}\) \(3678\)

[In]

int((d*x+c)^(1/2)*(f*x+e)^(1/2)*(h*x+g)^(1/2)/(b*x+a),x,method=_RETURNVERBOSE)

[Out]

((d*x+c)*(f*x+e)*(h*x+g))^(1/2)/(d*x+c)^(1/2)/(f*x+e)^(1/2)/(h*x+g)^(1/2)*(2/3/b*(d*f*h*x^3+c*f*h*x^2+d*e*h*x^
2+d*f*g*x^2+c*e*h*x+c*f*g*x+d*e*g*x+c*e*g)^(1/2)+2*((a^2*d*f*h-a*b*c*f*h-a*b*d*e*h-a*b*d*f*g+b^2*c*e*h+b^2*c*f
*g+b^2*d*e*g)/b^3-2/3/b*(1/2*c*e*h+1/2*c*f*g+1/2*d*e*g))*(g/h-e/f)*((x+g/h)/(g/h-e/f))^(1/2)*((x+c/d)/(-g/h+c/
d))^(1/2)*((x+e/f)/(-g/h+e/f))^(1/2)/(d*f*h*x^3+c*f*h*x^2+d*e*h*x^2+d*f*g*x^2+c*e*h*x+c*f*g*x+d*e*g*x+c*e*g)^(
1/2)*EllipticF(((x+g/h)/(g/h-e/f))^(1/2),((-g/h+e/f)/(-g/h+c/d))^(1/2))+2*(-1/b^2*(a*d*f*h-b*c*f*h-b*d*e*h-b*d
*f*g)-2/3/b*(c*f*h+d*e*h+d*f*g))*(g/h-e/f)*((x+g/h)/(g/h-e/f))^(1/2)*((x+c/d)/(-g/h+c/d))^(1/2)*((x+e/f)/(-g/h
+e/f))^(1/2)/(d*f*h*x^3+c*f*h*x^2+d*e*h*x^2+d*f*g*x^2+c*e*h*x+c*f*g*x+d*e*g*x+c*e*g)^(1/2)*((-g/h+c/d)*Ellipti
cE(((x+g/h)/(g/h-e/f))^(1/2),((-g/h+e/f)/(-g/h+c/d))^(1/2))-c/d*EllipticF(((x+g/h)/(g/h-e/f))^(1/2),((-g/h+e/f
)/(-g/h+c/d))^(1/2)))-2*(a^3*d*f*h-a^2*b*c*f*h-a^2*b*d*e*h-a^2*b*d*f*g+a*b^2*c*e*h+a*b^2*c*f*g+a*b^2*d*e*g-b^3
*c*e*g)/b^4*(g/h-e/f)*((x+g/h)/(g/h-e/f))^(1/2)*((x+c/d)/(-g/h+c/d))^(1/2)*((x+e/f)/(-g/h+e/f))^(1/2)/(d*f*h*x
^3+c*f*h*x^2+d*e*h*x^2+d*f*g*x^2+c*e*h*x+c*f*g*x+d*e*g*x+c*e*g)^(1/2)/(-g/h+a/b)*EllipticPi(((x+g/h)/(g/h-e/f)
)^(1/2),(-g/h+e/f)/(-g/h+a/b),((-g/h+e/f)/(-g/h+c/d))^(1/2)))

Fricas [F(-1)]

Timed out. \[ \int \frac {\sqrt {c+d x} \sqrt {e+f x} \sqrt {g+h x}}{a+b x} \, dx=\text {Timed out} \]

[In]

integrate((d*x+c)^(1/2)*(f*x+e)^(1/2)*(h*x+g)^(1/2)/(b*x+a),x, algorithm="fricas")

[Out]

Timed out

Sympy [F]

\[ \int \frac {\sqrt {c+d x} \sqrt {e+f x} \sqrt {g+h x}}{a+b x} \, dx=\int \frac {\sqrt {c + d x} \sqrt {e + f x} \sqrt {g + h x}}{a + b x}\, dx \]

[In]

integrate((d*x+c)**(1/2)*(f*x+e)**(1/2)*(h*x+g)**(1/2)/(b*x+a),x)

[Out]

Integral(sqrt(c + d*x)*sqrt(e + f*x)*sqrt(g + h*x)/(a + b*x), x)

Maxima [F]

\[ \int \frac {\sqrt {c+d x} \sqrt {e+f x} \sqrt {g+h x}}{a+b x} \, dx=\int { \frac {\sqrt {d x + c} \sqrt {f x + e} \sqrt {h x + g}}{b x + a} \,d x } \]

[In]

integrate((d*x+c)^(1/2)*(f*x+e)^(1/2)*(h*x+g)^(1/2)/(b*x+a),x, algorithm="maxima")

[Out]

integrate(sqrt(d*x + c)*sqrt(f*x + e)*sqrt(h*x + g)/(b*x + a), x)

Giac [F]

\[ \int \frac {\sqrt {c+d x} \sqrt {e+f x} \sqrt {g+h x}}{a+b x} \, dx=\int { \frac {\sqrt {d x + c} \sqrt {f x + e} \sqrt {h x + g}}{b x + a} \,d x } \]

[In]

integrate((d*x+c)^(1/2)*(f*x+e)^(1/2)*(h*x+g)^(1/2)/(b*x+a),x, algorithm="giac")

[Out]

integrate(sqrt(d*x + c)*sqrt(f*x + e)*sqrt(h*x + g)/(b*x + a), x)

Mupad [F(-1)]

Timed out. \[ \int \frac {\sqrt {c+d x} \sqrt {e+f x} \sqrt {g+h x}}{a+b x} \, dx=\int \frac {\sqrt {e+f\,x}\,\sqrt {g+h\,x}\,\sqrt {c+d\,x}}{a+b\,x} \,d x \]

[In]

int(((e + f*x)^(1/2)*(g + h*x)^(1/2)*(c + d*x)^(1/2))/(a + b*x),x)

[Out]

int(((e + f*x)^(1/2)*(g + h*x)^(1/2)*(c + d*x)^(1/2))/(a + b*x), x)